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Abstract 


We study the Boussinesq equation from the point of view of a multiple¬ 
time reductive perturbation method. As a consequence of the elimination 
of the secular producing terms through the use of the Korteweg~de Vries 
hierarchy, we show that the solitary-wave of the Boussinesq equation is a 
solitary-wave satisfying simultaneously all equations of the Korteweg-de Vries 

hierarchy, each one in an appropriate slow time variable. 
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I. INTRODUCTION 


As is well known, the Bonssinesq model eqnation 

Utt ^XX '^XXXX 3(11 ')xx 0 , ( 1 ) 

where u{x,t) is a one-dimensional held and the snbscripts denote partial differentiation, 
is completely integrable.0 It is considered as an intermediate long-wave equation since its 
long-wave limit with a further restriction to waves moving in only one direction yields, at 
the lowest order, the Korteweg-de Vries (KdV) equation 

Uf QuUx Uxxx 0 7 ( 2 ) 


which is an equation governing general weak nonlinear long-wave dynamics of dispersive 
systems-i Equation (|^) has N-soliton solutions. In particular, its solitary-wave solution is 
of the form! 


u = —2k‘^ sech^ k (x — Vl — 4:k^ t) 


( 3 ) 


where k is the wavenumber. Accordingly, its long-wave limit is related, also at the lowest 
order, to the solitary-wave solution of the KdV equation. 

In a recent work,@ we have considered a perturbative scheme based on the reductive 
perturbation method of Taniuti,i modified by the introduction of an infinite number of slow 
time-variables, which were given by rs = e^t, ts = e^t, etc. Then, we have shown that, as a 
consequence of a natural compatibility condition, a wave field satisfying the KdV equation in 
the time must also satisfy all equations of the KdV hierarchy,! each one in a different slow 
time variable. The main reason for introducing these time variables was that they allowed 
for the construction of a perturbative series, valid for weak nonlinear dispersive systems, 
which was free of solitary-wave related secularities. 

As stated above, the solitary-wave solution of the Bonssinesq equation tends to the KdV 
solitary-wave for small wave-numbers. What we will show here is that, by making use of 
the perturbative scheme with multiple slow time-scales, the solitary-wave of the Bonssinesq 
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equation may be written, in the slow variables, as a solitary-wave solution to the whole set 
of equations of the KdV hierarchy, each one in a different time-scale. This result follows 
both, from the general long-wave perturbation theory, and from the observation that the 
perturbative series truncates for a solitary-wave solution to the KdV hierarchy equations, 
rendering thus an exact solution for the Boussinesq equation. 

The paper is organized as follows. In Section II the multiple-time formalism is introduced 
for the Boussinesq equation, and the hrst few evolution equations are obtained. In Section 
III we discuss how the KdV hierarchy equations show up, and in Section IV we show how 
they can be used to eliminate the solitary-wave related secularities of the evolution equations 
for the higher-order terms of the wave held. In Section V, by returning from the slow to 
the laboratory coordinates, we obtain the above mentioned relation between the solitary- 
waves of the Boussinesq and the KdV hierarchy equations, namely that the solitary-wave of 
the Boussinesq equation may be written, in the slow variables, as a solitary-wave satisfying 
simultaneously all equations of the KdV hierarchy . And hnally, in Section VI, we summarize 
and discuss the results obtained. 


II. THE MULTIPLE TIME FORMALISM 


In order to study the long-wave limit of eq.(|^), we will introduce slow space and time 
variables based on the long-wave limit of the linear dispersion relation 

uj = k{l + ^ . (4) 

This limit corresponds to take 


k = en,, 


( 5 ) 


with e a small parameter. Expanding the dispersion relation (^), the solution of the corre¬ 
sponding linear Boussinesq equation is given simply by 


u = a exp i 




eK{x -t)- -^t + 


16 


-t + 


( 6 ) 
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Based on this solution, we define now the slow space coordinate 


^ = e{x-t), 


( 7 ) 


and the infinite sequence of slow time coordinates 


eH 


T3 = -- 


Accordingly, we have that 


T-S = 


eH 


Ty = - 


eH 

16 


( 8 ) 


and 


d d 
dx ^ 


d _ d e^ d e^ d e^ d 
dt ^ 2 Ots 8 dr^ IQdrj 


( 9 ) 


( 10 ) 


Notice that in the definition of r2n+i, we have already assumed specific slow time normaliza¬ 
tions, as inspired by the long-wave expansion of the dispersion relation. As we are going to 
see, these normalizations are exactly those necessary to cancel out the solitary-wave related 
secularities appearing in the higher order evolution equations. 

Returning to the nonlinear problem, we make now the expansion 


u = e^u = e^ [uq + e^U 2 + e'^u^ -!-■■■) , 


( 11 ) 


and we suppose that U2n = U2n{^, t's, Ts, ...), n = 0,1, 2 ,..., which corresponds to an extention 
in the sense of Sandri.l Substituting it, together with eqs.(|^) and (p! 0 [), into the Boussinesq 
equation (J^), the resulting expression, up to terms of order is: 

f \ e^ / 


_ 

4 I d^dr^ i 8 {d^dxY dr^dr^ 


2 , o,2„, , ,4 


U 


-3^ [(mo)^ + ‘2e^uoU2 + e'^{2uoU4 + (^ 2 )^) H-J = 0 . (12) 

We now proceed to an order-by-order analysis of the problem. At order e° we get 

= 0 . ( 13 ) 


A 
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Integrating once and assnming a vanishing integration constant, we obtain 


duo duQ d^uo „ 

- buo = 0 , 


9X3 




( 14 ) 


which is the KdV equation. 
At order e^, eq.(^) yields 


9 


du 2 9 , , duo 

■WF, - 


1 d^Un 1 9^Mn 


(15) 


4 9^9x5 4 9x3 2 

Using eq. (|l4|) , integrating once in ^ and assuming a vanishing integration constant, we obtain 

1 duQ 


du2 „ 9 9^M2 

^ - 6 ^(moM2 ) + ^ 


Id^uo d^uo 9duod'^uo 

+ 3mo Cl/-'? + TTWXXWTX)-9 (Mo) -ttt- ■ (16) 


/3 uq uq- ‘iui^ 4 9^3 ^^3 2 9 .^ 9.^2 

Equation ([T 6 |), as it stands, presents two problems. First, as Uqt^ is not known a priori, it 
cannot be solved for U 2 - In the next section we will show how to obtain the evolution of Uq 
in the time X 5 independently. The second problem is that the term (9 ®mo/ 9^^), as a source 
term for M 2 , is a secular producing term when uq is chosen to be a solitary-wave solution of 
the KdV equation. For instance, if we take the solution of eq. (0 proportional to [sech^ 9 ], 
then (9 ^Mo/ 9^^) will contain a term proportional to [sech^ 6 tanh^]. Being a solution of the 
homogeneous part of eq. ( 0 , this term produces a resonance, giving rise to non-uniformities 
in the perturbative series. It will turn out, however, that uqt^ can be used to cancel out this 
secular term. 


III. THE RISE OF THE KORTEWEG-DE VRIES HIERARCHY 

As we have seen, the held Mq satishes the KdV equation in the time X 3 : 

■^Ors = + 6 mo = F 3 . (17) 

The evolution of the same held Uq in any of the higher-order times X 2 n+i can then be obtained 
in the following way.S First, to have a well ordered perturbative scheme we impose that each 
one of the equations 
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M 0 t 2 „+i — F2n+l{Uo, Uq^, . . .) (18) 

be e-independent when passing from the slow {uq, ^,T 2 n+i) to the laboratory coordinates 
{u,x,t). This step selects all possible terms to appear in F 2 „+i(mo, .. .)• For instance, 

the evolution of uq in time rs is restricted to be of the form 


UQts — Ctrto(55) + + (/3 + 7 )Mo^MoC5 + , 


(19) 


where a, (3, 7 and 6 are constants. Then, by imposing the natural (in the multiple time 
formalism) compatibility condition! 


('^Oraj (^0r2n+l) : 

^ '' T2n+1 ^ ^ / ra 


( 20 ) 


or equivalently. 


(-^3)t-2„+i - i.^ 2 n+l)rs > (21) 

with Fs given by eq. (|T7|) , it is possible to determine any F 2 n+i- As it can be verified,! the 
resulting equations are those given by the KdV hierarchy. In particular, for mqt-s and 
we obtain respectively 

wot5 = Uo{50 ~ — 20'Uog'Wog^ + 30 (mo)^mo$ , (22) 


and 


U0t7 — — rto{70 + 14 moMo{5$) + 42mo^Mo(4^) + 

+ — 280«oMo^mo^^ — 70(mo^)^ — 70(mo)^mo5cc 


(23) 


In principle, one could have an arbitrary constant multiplying the right-hand side of eqs.(p^ 
and (p3D, which would correspond to an arbitrarity in the slow time normalizations. However, 
as we will see in the next section, the dehnition of the slow time variables we took implies 
that these constants must be chosen to be one, since in this case the perturbation theory 
is automatically rendered free of secularities. This choice is also the one which makes the 
linear limit of the perturbation theory compatible with the linear theory coming directly 
from eq.(|]). 
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IV. HIGHER ORDER EVOLUTION EQUATIONS 


From this point on, we are going to consider some specific solutions to our equations. 
First of all, we assume the solution of the KdV equation (|T7|) to be the solitary-wave solution 


Uq = —2K^sech^ + 9 


(24) 


where 6^ is a phase. Moreover, since mq must satisfy also the equations of the KdV hierarchy, 
we assume that uo given by be also a solitary-wave solution to all equations of the KdV 
hierarchy, each one in a different slow-time variable. This means that Uq is actually 


Mo = —2K^sech^ K,^ — + 


(26) 


We return now to eq. ([l6|) for U 2 - Substituting mqtb from eq. (|^) , we obtain 


U2 t 3 — Q{UoU2)^ + — - ~3(Mo)^'Wo^ + ~ 


(26) 


In passing we notice that the use of the KdV hierarchy equation to express mqts automatically 
canceled out the secular-producing term Mo(50 - Moreover, using the solitary-wave solution 


for Mo, we see that the right-hand side of eq. (|^) vanishes, leading to 


U2t3 — 6(MoM2)^ + U2^^^ — 0 , 


(27) 


which is the linearized KdV equation. We will assume for it the trivial solution 


M 2 = 0 . 


With this result, order is solved for the particular case we chose. 
At order e^, and already assuming that M 2 = 0, eq.(|^ gives 


1 r 


M4t3^ — 6(moM 4)^^ + M4(4g) — - —Mot^^ + Mq 


JT 3 T 5 


(28) 


(29) 


Using equations (p!?]) and (p2|) to express Morg and Morg respectively, and integrating once in 
we obtain 
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M4t3 ~ 6(moM4)5 + — o ~ W 077 — '^0(70 + 16moMo( 5C) ~ 90('*^o)^'*^occc + 

o '- 

+ 70Mo^^Mo5gg + 40MogMo(40 “ SOQuoUq^Uq^^ + 180(mo)^Moc - 60(Mog)^] . (30) 

The term mo( 75 ) is the only resonant, that is, secular producing term to the solution U 4 . 
Then, in the very same way we did before, we use the KdV hierarchy equation (^) to 
express Kqty- After we do that, the secular producing term is automatically canceled out. 


and eq.f^OD becomes 


M4t3 — QiuoUi)^ + — — 2[mo^Mo( 4^) — Mo'Wo(5|) + 

+ - 5(mo^)^ + 10(mo)^Mo 555 - 20(Mo)^Mod • (31) 

Substituting again the solitary-wave solution (^) for uq, we can easily see that the nonho- 
mogeneous term of eq. ( 0 ) vanishes, leading to 

M 4 T -3 — 6('160^4)^ + ^ 4 ??? ~ 9 • ( 92 ) 


And again, we take the trivial solution 


M 4 = 0 . (33) 

It is easy to see that this is a general result that will repeat at any higher order: for n > 1, 
the evolution of U 2 n in the time Ts, after using the KdV hierarchy equation to express Mor 2 „+i 
and then substituting the solitary-wave solution (^5]) for uq, is given by a homogeneous 
linearized KdV equation. Consequently, the solution 

U2n = 0 ; n > 1 (34) 


can be assumed for any higher order. 


V. BACK TO THE LABORATORY COORDINATES 

Let us now take the solitary-wave solution to all equations of the KdV hierarchy. 
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Uq = —2K^sech"^ — 4:K'^T3 + — 6Ak tj + 


( 36 ) 


and rewrite it in the laboratory coordinates. First, recall that we have made the expansion 


u = e^u = [uq + e^U 2 + e'^M 4 + ■ ■ ■) 


(36) 


Thereafter, we have fonnd a particnlar solntion in which 


M2n = 0 ; n > 1. 


(37) 


Conseqnently, expansion (^6l) trnncates leading to an exact solntion of the form 


U = e Uq, 


(38) 


with Uq given by eq.(^). Moreover, from eq.(^ we see that the wave-nnmber k is written 
in terms of the corresponding laboratory one by 


K = e . 


(39) 


Finally, the slow coordinates ^ and T 2 n+i are related to the laboratory ones, x and t, according 
to eqs. (0) and (H). Then, in the laboratory coordinates, the exact solntion (|38D is written as 


u = —2fc"^sech"^ k 


X 


(l - 2P - 2/t* - 4/t“ + • • ■) t 


(40) 


Now, the series appearing inside the parenthesis can be snmmed, with the resnlt 


l-2k‘^ -2k^-Ak^+ ■ ■■=(!- Ak'^) 


1/2 


(41) 


Therefore, we get 


u = —2/c^sech^ 


k[x — \/l — 


(42) 


which is the well known solitary-wave solntion of the Bonssinesq eqnation (|^). 
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VI. FINAL COMMENTS 


By applying a multiple time version of the reductive perturbation method of Taniutill 
to the Boussinesq model equation, and by eliminating the solitary-wave related secular 
producing terms through the use of the KdV hierarchy equations, 0 we have suceeded in 
establishing a relation between the solitary-wave satisfying all the equations of the KdV 
hierarchy and that of the Boussinesq equation. In other words, we have shown that the 
solitary-wave of the Boussinesq equation is given, in slow variables, by the solitary-wave 
satisfying simultaneously all the equations of the KdV hierarchy. Accordingly, while the 
KdV solitary-wave depends only on one slow time variable, namely rs, the solitary-wave of 
the Boussinesq equation can be thought as depending on the inhnite sequence of slow time 
variables. 

The above considerations put in evidence the universal character played not only by 
KdV, but by all equations of the KdV hierarchy in relation to general weak nonlinear dis¬ 
persive systems. For such systems, as we have already said, it is always possible to dehne 
slow variables in which the KdV equation emerges at the lowest relevant order of the reduc¬ 
tive perturbation method, and consequently, compatibility and secularity-free requirements 
imply that all equations of the KdV hierarchy emerge as well. In the case where the pertur¬ 
bative series truncates, we may then obtain an exact solution of the original equation, which 
is, in this sense, reconstructed from the perturbative expansion. Of course, the Boussinesq 
solitary-wave is a well-known solution, but not the least, it remains that a possible method 
of construction of solutions for more involved system of equations can be envisaged. The 
return to the laboratory coordinates, then, make the connection between a solution of the 
KdV hierarchy and that of the equations governing the original system. 

To conclude, we may conjecture that, whenever the original model equation has an 
exact solitary-wave solution, the series may somehow be truncated (or eventually summed) 
and a relation will be established between the solitary-wave of the KdV hierarchy and 
that of the original equation. On the other hand, when the original nonlinear dispersive 


10 


system does not present an exact solitary-wave solution, the series will not truncate. In 
this case, a secular-free expansion can still be obtained and the process of returning to 
the laboratory coordinates can be made order-by-order at any higher order, implying in a 
sucessive solitary-wave velocity renormalization.&i 
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